One Dimensional Continuum 
Falicov-Kimball Model in the Strongly 

Correlated Limit 



■r::j- ! Robert James Bursill 



Department of Mathematics 

The University of Melbourne 
Emaih r.bursill@sheffield.ac.uk 

Short Title - Continuum Falicov-Kimball Model 
Classification Number - 05.30 - Statistical Mechanics — Quantum 



(N 



> 

in 
vn 

o 

1^- [ Keywords - quantum statistical mechanics, model system, interacting fermion model 

a 

B 

a 
o 
o 

> 



Abstract 

In this paper we study the thermodynamics of the one dimensional continuum analogue 
of the Falicov-Kimball model in the strongly correlated limit using a method developed by 
Salsburg, Zwanzig and Kirkwood for the Takahashi gas. In the ground state it is found 
that the / electrons form a cluster. The effect of including a Takahashi repulsion between 
/ particles is also studied where it is found that as the repulsion is increased the ground 
state / electron configuration changes discontinuously from the clustered configuration to 
a homogeneous or equal spaced configuration analogous to the checkerboard configuration 
which arises in the lattice Falicov-Kimball model. 



1 Introduction 

The Falicov-Kimball model m was initially proposed as a simple model for the semiconductor- 
metal transition in SmBg and transition metal oxides. The model consists of two species of 
electrons — localised / electrons and a itinerant band of d electrons. The Hamiltonian is 

n = Y,ekclcu + UY,d\d,flU (1-1) 

k i 

where fi is the destruction operator for / electrons on site i, di is the destruction operator for 
the d-band Wannier state at site i and Ck is the destruction operator for the d electron Bloch 
state with wave vector k and kinetic energy Cfe. 

The only interaction then is an on-site interaction between / and d electrons. Of course 
the model could be extended to include longer range f-d interactions, /-/ (lattice gas type) 
interactions and d-d interactions. Spin could be introduced for the d electrons with Hubbard-type 
interactions. In [EJ and [pj the model was independently introduced as a model for crystallization 
and also viewed as a Hubbard model H where one of the spin species of electrons is immobile. 

A great deal more is known about the Falicov-Kimball model than is known about the Hub- 
bard model. In |ll| the model was treated within a mean field framework. It was deduced that for 
large enough U there is a first order phase transition which was interpreted as a semiconductor- 
metal transition. 

In p] it was established by means of operator inequalities that for the bipartite lattices with 
equal numbers of / and d electrons and a half-filled band, the ground state is ordered in any 
dimension, the / electrons forming a checkerboard pattern. For d > 2 it was shown further that 
the crystalline order persists at finite temperature. 

In (5| and ||6| sharp upper and lower bounds are obtained for the ground state energy in 
the two dimensional case. For some special values of the parameters the bounds are shown to 
coalesce and so the exact ground state energy is derived. Some of the results of |g] and |^ are also 
independently reproduced. In and g the thermodynamics of the model is investigated for 
large dimensionality d. An exact solution is obtained in the d = oo limit for the thermodynamic 
potential, the critical temperature, the order parameter and several correlation functions. 

In B the ground state of the one dimensional model is studied by means of second order 
perturbation theory and exact numerical calculations for periodic configurations of the / elec- 
trons. A restricted ground state phase diagram with a fractal structure is derived. In pO| the one 
dimensional model is studied for large negative U, the / electrons being interpreted as nuclei. 
An exact formula for the leading order contribution (for large \U\) to the ground state energy is 
derived from which it is deduced that the system forms atoms {d-f pairs) which have an effective 
repulsion between them if \U\ is sufficiently large. 

In this paper we study the one dimensional continuum analogue of the Falicov-Kimball model 
in the strongly correlated limit. Although in terms of mathematical complexity this model is a 
long way removed from the lattice model with finite interaction strength in higher dimensions, 
there may be some generally applicable physics which can be deduced from this study. We 
obtain the exact thermodynamic potential via a method first used by Salsburg, Zwanzig and 
Kirkwood fO] for the Takahashi gas ||l2|. We discuss the ground state energy as well as a 
physical interpretation of the model in terms of the Takahashi gas. Finally, we discuss the effect 
of introducing a Takahashi type interaction between the / particles. 

2 Continuum Model and Exact Solution in One Dimension 
in the Limit of Strong Correlations 

In the continuum analogue of the Falicov-Kimball model, we have a system of fermions {d elec- 
trons) and stationary scatterers (/ electrons). If there are M scatterers at positions yi, ■ ■ ■ ,yM 



then the equation for the d particle eigenfunction ip{x) is 



1 ^' 



i=i 



ijjix) = EiPix) 



(2.1) 



where E is the energy eigenvalue and we have scaled the fcrmion mass and Planck's constant 
appropriately. 

In the strongly correlated limit {U = oo) the one dimensional continuum analogue of the 
Falicov-Kimball model reduces to a system containing a number of impenetrable barriers con- 
figured along the line. That is, if the length of the line is V and there are M / electrons at 
positions < yi < y2 < ■ • • < Vm < V then the d electrons experience impenetrable barriers 
at positions yo = 0, j/i, . . . , j/m, and i/m+i = V|j or infinitely high wells in the regions (j/o, yi), 

{yi,y2),---,{yM,yM+i)- 

In this case the eigenvalue problem ( |2.l| ) is trivial, reducing to 

1 



— V'"(a;) = E^{x) 



(2.2) 

^(%) = for j=0,l,...,M + l (2.3) 

i;{x) ^ for X < and X > y (2.4) 

The single particle wavefunction for a state with wave number k > 1 m well 1 < j < M + I is 



fpjkix) 



,/^^sin (^hiHizVi^) y <x<y, 

otherwise 







and the corresponding energy eigenvalue is 



fc2 



(yj -j/i-i) 



(2.5) 



(2.6) 



To calculate the thermodynamic potential for this model, we begin by defining a canonical 
partition function for the system with M f electrons and N d electrons viz 



M+l 00 



Zmn = / ^y H H X! exp{-f3njkEjk) 
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(2.7) 



j=l fc=lrijfc=Oa 

where the primed sum is restricted to configurations with N d particles ie 

M+l 00 

j=i fc=i 



and 



Vy=\ dyu / dyu-i ■ ■ 
Jo 
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dyi 



(2.8) 



(2.9) 



The restriction on the sum can be removed by defining a grand canonical partition function 
Q{V) = ^z^Zm^v (2.10) 



N=l 



/M+l 00 
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(2.11) 



^We assume fixed end boundary conditions. 



where use has been made of (2/7) 
Making the definition 
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f + ze-l''^"!-" 
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we have, using (|2!|), (|j|), ( ^.ffD and ( ^.f2D that 



„ M+l 



(2.12) 



(2.13) 
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nV2-Vi)nVi)dyi (2.f4) 
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(2.15) 



where /*■■■*/ denotes the (Af + f )-fold Laplace convolution of /. 
We take the Laplace transform of Q with respect to V viz 



/>oo 

Cq{s) = / e-'^Q{V)dV 
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'e-^^(/*...*/)(F) 



[F{s)] 



M+l 



where 



F{s) = / f{x)e-'''dx 



(2.16) 

(2.17) 
(2.18) 

(2.19) 



denotes the Laplace transform of / and we have made use of ( 2.15| ) and the Laplace convolution 
theorem. 



Applying the Laplace inversion formula to ( 2.1g ) we then have 



-1 /'C^lOD 

Q{V) = ;— / [F{s)] 
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(2.20) 



where the contour lies to the right of all singularities of F{s). 

The singularity structure of F can be deduced from a different representation of / which we 
derive in appendix H viz 
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where 



Xi 



poo 

/ log n + ze-l^^'j dx 



denotes the thermodynamic potential for a system of free d particles, 
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(2.21) 
(2.22) 

(2.23) 
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The Uk and Vk form increasing sequences so F{s) has simple poles at 



s = Xi 

s = Xf - 27r (wfc ± TOfc) fc = 1,2,3,... 

s = Xf- 47rwfe A: = 1,2,3, . . . 

s = Xf — 27r (ufe ± iuk + vi ztiui) 1 < k < I 



(2.25) 



In ( 2.20| ) then, we must take c > Xf- 

Now the thermodynamic potential x is given by 



lim 



log Q{V) 

V 



where in taking the limit, the well density 
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(2.26) 



(2.27) 



is held constant. We evaluate the right hand side of ( 2.20 ) for large V and M by the method of 
steepest descents. We write 
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exp [VH{s)] ds 



c> Xi 



(2.28) 



where 



H{s) = s + mlogF{s) (2.29) 

and so, assuming that H{s) has a single stationary point sq in the half plane Sfts > Xf we have 

X^Hiso) (2.30) 



From the pole structure of F (2.25) and the fact that x is real, we deduce that sq lies on the 
real axis, the schematic form of H{s) being given as in Fig. |l|. That is, H{s) has a branch point 
at s = Xf ^'^d H{s) ^ s as s ~> oo. sq is therefore a local minimum of H{s) on the real axis viz 



X — min H{s) 

s>X{ 



(2.31) 



The steepest descents method was first employed by Salsburg, Zwanzig and Kirkwood u^ 
in their calculation of equilibrium properties of the Takahashi gas [|l2[ which, as we shall in 
section 0, has some similarities with the present model. 

The thermodynamics of the model can be determined from x- For instance, the fugacity can 
be rendered as a function of the temperature, the d electron density n = N/V and the / electron 
(well) density m by inverting the relation 



dx 



The Helmholtz free energy is then given by 



V' = M ^ 

np 



(2.32) 



(2.33) 



In general, sq cannot be calculated analytically. In the limit /3 — > oo however, the function H{s) 
simplifies and the ground state energy can be derived. 



3 Ground State Energy 

The ground state energy (per d electron) is given by 



where 



PG 

Eg = A* 

n 



PG = lim -s 



(3.1) 



(3.2) 



denotes the ground state pressure and in (|3.1| ) fi is determined as a function of the d electron 
density by inverting the limiting form of (2.32|) viz 



dpG 
9/i 



To find pq we require the limiting form of H{s)/ j3. Now, using (2.22) we have 



Xf - r\og{l + e^i^--') 



dx 



^ j3 I {fJ- ~ x'^) dx 
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Writing s = Pfi^^^t then, we use ( |2.19| ) and (|t|) to get 

F{s) = / TT { 1 + cxp [l3fi (l - k'^/fix'^)] } exp [~/3fit^x] dx 
Jo ,._, 
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1 f 

= -— TT { 1 + exp [/3/i (1 - fc Va;^)] } exp [^Pfitx] dx 



fe=i 
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where 
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and we have used the identity 
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Clearly g{x,t) ~ (2/3 — t)x as x ^ oo and hence the integral ( B.IOD converges if i > 2/3 as we 
would expect fr om (|3.6| ) and the requirement that s > Xf- 
Evaluating (B.IOD by Laplace's method we then have 



F{s) ~ exp 



l3^Ta.axg{x,t) 

x>0 



whence from ( 2.29 ), ( 2.31 ) and (3^) we obtain 



Pg — min hit) 

t>2/3 



(3.13) 



(3.14) 



where 



h{t) ^ lim ^^; ' (3.15) 

= fi ' t + mijLinaxg{x,t) (3.16) 

a;>0 

g{x, t) is sketched as a function of x for various values of t in figureH. We see that maxg(a::, t) = 

x>0 



g{0,t) — for aU t > 2/3. From ( 3.16 ) we therefore have 

h{t) = fi^^^t (3.17) 



and so from (3.14) we have 



Pg = |m^/^= limf EEpg) (3.18) 

That is, as /3 — > oo, So ~ Xf + 0(1), the large /? form of H{s) being as indicated in Fig. H. 

Pg therefore equates to Pq — the ground state pressure for a free system of d electrons. Using 



(3.3) we get 

li = n^ (3.19) 

and hence the ground state energy is 

SG = £g (3-20) 

where Eq denotes the well known free fermion ground state energy 



.(f) "■ 



(3.21) 



The ground state is therefore a state where all of the / particles are clustered at the ends and 
all but one of the wells is unoccupied and of zero width, the d electrons occupying the large well. 
In the following section we interpret this result by considering the analogy between the present 
model and the Takahashi gas. 

4 Physical Interpretation and Inclusion of a Takahashi In- 
teraction Between / Electrons 

In order to physically interpret the above result we recall that, for a Takahashi gas Q, we have 
a system of M particles at positions < j/i < . . . < j/j\/ < V with a repulsive interaction 

J2 <t>{yj-yj-i) (4.1) 

0<i<j<M+l 

where 

oo < X < a 

(f){x) = { v{x) a < X <2a (4.2) 

x> 2a 



and in (4.2) a > is the hard core exclusion distance and v{x) > 0. 

Because of the hard core exclusion, the interaction only acts between nearest neighbours and 
the partition function is 

Jq Jo 

y2 

g-/30(a2-ai)g-/30fei)^y^ (4.3) 





As mentioned in section g, Laplace convolution and the method of steepest descents have been 
employed |11[ t o solve for the Helmholtz free energy of the Takahashi gas. 



Taking (4.3) as a starting point, we could define a model partition function with longer range 
forces by relaxing the requirement that (j>{x) vanishes for x > 2o. As pointed out however |K| 
such a partition function is unphysical in that only nearest neighbours interact and the specific 



two-body potential form (4.2) is required for (4.3) to be valid. 

We note however that the grand canonical partition function from the present model (2.14) 

That is, the nearest neighbour potential 



(x) with i^g/M. 



is of the form (O) if we identify 

between / ele ctron s can be considered an effective interaction induced by the d electrons. As can 
be seen from ( 2.12 ), this interaction is temperature dependent, long ranged and attractive. This 
sheds some light on the clustering of the / electrons observed in section (^) . A further physical 
interpretation arising from an approximate evaluation of the thermodynamic potential is given 
in appendix ^. 

Intuitively, from the results mentioned in section |lj for the lattice Falicov-Kimball model we 
might expect that the ground state / electron configuration would be the analogue of the checker- 
board configuration with the / particles homogeneously spaced along the line. In appendix y 
we derive the ground state energy from the equally spaced configuration 
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(4.4) 



We show explicitly that this lies above the exact ground state energy (3.21) which arises from 
the highly inhomogeneous / electron configuration mentioned in sectionM 

It is natural then to investigate the effect of including a Takahashi term (E|2) (a short range 
repulsive potential between / particles) in the model. In this case the barriers cannot cluster 
together and the resulting model should more closely resemble the lattice model where the / 
electrons cannot cluster into a very small region. 



With the Takahashi term included, the grand canonical partition function has the form ( 2.14 ) 
but with / now defined by 



/(-) - n [ 

k=l 
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^I3k^/x' 



-04,(x) 



(4.5) 



The pole structure of F{s) does not change and the thermodynamic potential is still obtained 
from (|3^). 

We restrict our attention to the Tonks case jlj] where v{x) = in order to calculate the effect 
that a repulsive /-/ interaction has on the ground state. As the hard core exclusion distance a 
is increased from to the extreme value 1/m, the number of allowed barrier configurations is 
reduced until the only allowed configuration is the homogeneous configuration. The ground state 
energy sq should in turn increase monotonically from the free fermion value (3.21) to the equal 
spacing value (4.4). Alternatively, the ground state pressure pq should decrease monotonically 
from the free fermion value (3.18) to the equal spacing value 
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(4.6) 



which is derived in appendix |C[ 



To derive the ground state pressure we repeat steps (3.7)-(3.10) with the modified form (4.5) 
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F{s) = / n{l + exp[/3^(l-A:V^a;2)]}exp[-/3/it77Ia;]e-'^*(^)da; (4.7) 

■^0 fc=l 

= / Y[{l + cxp[f3fi{l~k^/fix'^)]}exp[-Pnt^x]dx (4.8) 



-— Jl {l + exp [Pfi (l - k'^/x'^)] } exp [-/?/ite] dx 



Applying Laplace's method once more we have 



Pfi max g{x,i) 



F{s) ~ exp 
and we again have ( 3.14 ) but with ( |3.16| ) generalised to 



h{t) ~ fi' t + mil max g{x,t) 

x>ay^ 



(4.9) 
(4.10) 

(4.11) 
(4.12) 



Clearly h(t) and therefore pc is a monotonically decreasing function of a. As a approaches 
1/m we find that 



max g{x,to) 



g{ay/Jl,to) 
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(4.13) 
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1|_V^ (4.15) 



where tg 



lim so/f3fi^''^ is the minimizing value in (3.14). As expected, it follows from (4.6) 
(3— ►oo ^ ' ^ — ' 



and (3.14) then that pQ 



P^^ 



1/m. 



Numerical investigations of pc as a function of a reveal that the decrease from Pq to Pq 
is not smooth. We illustrate this for the case where fj, = 5 and ?7i = 2 in Fig. g where we plot 
h{t) versus t for various values of a and Fig. H where we plot where we plot pq as a function of 



(ES) 



We 



a. We see that there is a critical value of the exclusion Oc above which pQ equates to p^ 

are naturally lead to the following characterization of the ground state / electron configuration. 

For a — the / particles lie in the highly inhomogeneous configuration where all wells 
have zero width except one which has width V. For < a < Oc all wells must be at least 
as wide as a but inhomogeneity in the / configuration persists with one well wider than the 
others. For Oc < a < 1/to, the ground state configuration is the homogeneous or equal spaced 
configuration, each well having width 1/m, even though inhomogeneous configurations are allowed 
for flc < a < 1/to,. 

We note from Fig. H that pq may in general have further singularities in the region < a < Ac. 
These singularities relate not to the well configuration but to the distribution (or filling) of the 
d electrons in the wells. For sufficiently small a, the large well has width V — Ma (all the small 
wells have the minimal width a) and all of the d electrons lie in the large well. That is, the 
ground state pressure is that of a free fermion system with a reduced volume. It is easy to show 
that in this case the ground state pressure is (1 — ma)pQ . As a is increased towards a^, critical 
values occur at which it becomes energetically favourable to fill the smaller wells. In the specific 
case depicted, there is a single critical value of the exclusion ac € (0, ac) such that pQ equates to 
(1 — ma)pQ for < a < Oc and filling of the smaller wells occurs for Oc < a < ac. 



5 Discussion 

In this chapter we have studied the thermodynamics of a one dimensional continuum analogue of 
the Falicov-Kimball model in the limit of strong correlation. As mentioned in the introduction, in 
terms of mathematical and physical complexity, this model is a long way removed from a lattice 



model with finite strength interactions in higher dimensions . As shown in section ^ however, 
the exact solution does afford a neat physical interpretation and the addition of a repulsive /-/ 
interaction enriches the model leading to a ground state analogous to that observed in lattice 
models if the repulsion is sufficiently strong. Finally, we indicate some directions for further 
work. 

The calculation of the ground state / electron configuration and d electron distribution in 
the case where the /-/ interaction is of the Tonks form could be taken further and made more 
precise. Also, a study of the ground state when v{x) is non-zero may prove interesting because 
in this case the Takahashi repulsion will contribute directly to the ground state energy and not 
simply through the reduction of allowed / configurations. It should be possible to generalise 
to the present model a method used by Salsburg, Zwanzig and Kirkwood in the calculation of 
correlation functions for the Takahashi gas [O. The calculation of the /-/ correlation function, 
especially at zero temperature may be of some interest. It may be of interest to generalise the 
model to higher dimensions d where the barriers are impenetrable d— 1 dimensional hyperplanes. 
Finally, it might be interesting to investigate this model in the case where the d particles have a 
different energy spectrum, different statistics, or both. 

The author was supported by an Australian Postgraduate Research Award throughout the 
course of this work. The author acknowledges the encouragement of his Ph. D. supervisor, 
Professor Colin Thompson and is thankful to Dr. Bill Wood for some useful discussions. 

Appendix 

A A Representation for / 

In this appendix we derive the representation (|2.21) for the function /. Now from (2.12) we have 



oo ^ 

log/(x) = - Yl Hk/x)~-log{l + z) (A.l) 
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where 



T{w) = log (l + 2e~^"'') (A.2) 

J^{w) is absolutely integrable so the Poisson summation formula can be applied to ( |A.l| ) viz 

log/(x) = I E •^(2^M-^log(l + ^) (A.3) 

fc — — OO 

oo . 

= xxi + xYT{2TTkx)--\og{l + z) (A.4) 

where 

/"OO 

T{y) = / T{w)e''"ydw (A.5) 



denotes the Fourier transform of T and we have noted using (2.22) that T{0) — 2xi- 
On integrating by parts we find that 

T{y) = / T' {w)e''"y dw (A.6) 

Now 



so, above the real axis —T'{w) has simple poles at w — Wi = ±ui + ivi I = 1,2,3, .. . where ui 
and vi are the positive, real solutions of 



P{w^f = /3/i ± (2Z - l)7ri 



That is, the intersection of the two hyperboli 
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(A.8) 

(A.9) 
(A.IO) 



The explicit values of ui and vi are given by ( |2.23D and (2.24). The residue of —T'(w)e™y at 
w = Wi is — e™i y . Assuming y > 0, using ( [A.6| ) and closing the contour above the real axisQwe 
therefore have 

T{y) = --^ Y^ U'^tv + e'^^^A (A.ll) 
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Combining (A.4) with (A.ll) we arrive at 

log fix) = xxf - 2 log(l + ^) - E ^ E {e'^''"''^ 

fc=i 1=1 

Using the Taylor expansion 

ET = -i°g(i-'^) 

we then have 



^27Tikw, X 
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log fix) = xxf - - log(l + z) + E log (l - e^-™^) ( 
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(A.12) 



(A.13) 



(A.14) 



from which the required representation (2.21) follows 



B Approximate Low Temperature Calculation of the Ther- 
modynamic Potential and Further Physical Interpreta- 
tion 

In this appendix we give an approximate direct evaluation of the thermodynamic potential which 
should be reasonable at very low temperatures and which admits a neat physical interpretation. 
Now using ( 2.21 ) and ( 2.19 ) we can write 
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where 
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dx 



(B.l) 



(B.2) 



is analytic for ^s > 2ttvi. 

The right hand side of ( 2. 20] ) is dominated by the pole at s = Xf for large M. The terms 
arising from the other singularities of F are exponentially small. Using (B.l) we write 
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(B.3) 



^It is easily seen that tlie integral on the closing contour approaches zero as long as the closing contour is 
chosen to bisect the poles so that the denominator of —J^'{w) is bounded away from 0. 
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From ( 2.24 ) and ( 2.25| ) we see that the poles of F{s) accumulate on the line 5fts — xf as 
/3 — > oo. We therefore make the approximation G{s) ~ G (xt) and so from (B.3) we approximate 
the residue of {F{s))^'^~^^e'^^ at s = Xf by 
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V M+l 



(1 + ^) 2 j-^i 



M + l\[G{xi)f+'''V^- 



In this approximation then the dominant contribution to the partition function is 

2(1/) « ^^ 

Making use of the identity 
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where we have used the standard limit 
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Combining (B.5) and (B.IO) we then have 
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(B.4) 



{B.5) 



(B.6) 



(B.7) 
(B.8) 



(B.9) 
(B.IO) 

(B.ll) 
(B.12) 



Taking the thermodynamic limit (1/ ^ oo with m = [M + \)/V fixed) we obtain the approximate 
thermodynamic potential 



X ~ TO — m log TO + Xf + TO, log 



1 + mG (xf) 



(B.13) 



The above expression has a simple interpretation. The first term is an entropy term arising 
from the barrier configurations. The second term is the free fermion term arising from states 
where all wells have zero width except one with width V . The last term arises from states where 
a macroscopic number of wells have non-zero width. 
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Finally, from (B.2) and (2.21) we have 



G(xf) 
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^2iT{-Vk+iUk)x 



dx 



(B.14) 



Using methods from section || it can be shown that the la st ter m in ( B.13 ) does not contribute 
to the ground state pressure and so we recover the result (3.20) within this approximation. 



C Ground State Energy and Pressure for the Homoge- 
neous / Electron Configuration 

In this appendix we calculate some ground state quantities for the homogeneous or equal spaced 
/ configuration. That is, the position of the jth barrier is 



V 



Vj =J 



M + 1 m 



= — J = 1,- 



,M 



(C.l) 



In this case the M -\- 1 wells all have width \/m and so the energies (2^) reduce to 

Ejk = rr?k^ (C.2) 

For N d electrons, the lowest energy is achieved by filling the lowest energy states of the wells 
with as close to equal filling as possible. We write 



where 



is the integral part of n/m and 
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(C.3) 

(C.4) 

(C.5) 



The first q levels are filled in each well and the {q + l)th level is filled in {M + l)r of the wells. 
Using (C.2) then, the ground state energy per d particle is 
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(C.6) 



(C.7) 



where we have made use of ( 3.12| ). Using (C.4), ( |C.5| ) and a little manipulation, we then arrive 
at (O). 



We next show explicitly that the exact ground state energy (3.20) (from the highly inhomo- 
geneous configuration) lies below the ground state energy from the homogeneous configuration. 
Using ( pO| ), ( pH ), (g^), (^ and (|C||) we have 
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(C.8) 

(C.9) 

(CIO) 
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(C.12) 
(C.13) 
(C.14) 

(C.15) 

(C.16) 
(C.17) 



as required. 

Finally, we calculate the ground state pressure from the homogeneous configuration in terms 
of m and /i. To do so we find the grand canonical partition function for the specific barrier 
configuration (C^). The canonical partition function is 
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(C.18) 



and, using ( |C.2| ), the associated grand canonical partition function is 
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The thermodynamic potential is 



^(ES)(y) = lini 



log S(ES) 



V^oo V 



^ ^ log { 1 + exp [(3fi (1 - m'^k'^/fi)] } 



fc=i 

k=l 



mf3fi 






(for large (3) 
^ '1/2 



VT* 



(C.19) 
(C.20) 

(C.21) 

(C.22) 
(C.23) 

(C.24) 
(C.25) 



so the ground state pressure Pq — lim ^-g— is given by (4.6) 



,(ES ) 
/3^-^oo /^ 
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Figure 1: The schematic form of H(s) on the real axis (soUd Une) and the large /3 form (dot-dashed line). 



Figure 2: g(x, t) as a function of x for various values of t. 



Figure 3: h{t) as a function of t for various values of the exclusion distance a in the case where /i = 5 and 

m = 2. 



Figure 4: Ground state pressure pq (solid line) as a function of the exclusion distance a in the case where fi = 5 

(t) 
and m = 2. pQ equates to (1 — ma)p)^, (dot-dashed line)— 

with restricted volume V — Ma — in the region < a < dc- 



If) 
and m = 2. pQ equates to (1 — ma)p)^, (dot-dashed line) — the ground state pressure from a free fermion system 
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